On an integrable two-component Camassa-Holm shallow water system 
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The interest in the Camassa-Holm equation inspired the search for various generalizations of 
this equation with interesting properties and applications. In this letter we deal with such a two- 
component integrable system of coupled equations. First we derive the system in the context of 
shallow water theory. Then we show that while small initial data develop into global solutions, for 
some initial data wave breaking occurs. We also discuss the solitary wave solutions. Finally, we 
present an explicit construction for the peakon solutions in the short wave limit of system. 



PACS numbers: 05.45.Yv, 47.35.Bb, 47.35.Fg, 47.35.Jk 

In recent years the Camassa-Holm (CH) equation [l[ 

Wit + w«i + 2mu x + m x u — 0, m = u — u XXl (1) 

(u> being an arbitrary constant) has caught a great deal 
of attention. It is a nonlinear dispersive wave equation 
that models the propagation of unidirectional irrotational 
shallow water waves over a flat bed [l|, |2, lla, IJ] > as well 
as water waves moving over an underlying shear flow H . 
The CH equation also arises in the study of a certain 
non-Newtonian fluids [6( and also models finite length, 
small amplitude radial deformation waves in cylindri- 
cal hyperelastic rods Ml. The CH equation has a bi- 
Hamiltonian structure [8( (and an infinite number of con- 
servation laws) , it is completely integrable (see [l[ for the 
Lax pair formulation and 0] for the direct /inverse scat- 
tering approach), and its solitary wave solutions are soli- 
tons [H M, EI Gl IH with stable profiles [II El]. The 
equation attracted a lot of attention in recent years due to 
two remarkable features. The first is the presence of so- 
lutions in the form of peaked solitary waves or 'peakons' 



equation has smooth solutions which develop singulari- 
ties in finite time in the form of breaking waves 
- the solution remains bounded while its slope becomes 
unbounded in finite time [l9| ) as well as that of smooth 
solutions defined for all times [151 ] . These two phenomena 
have always fascinated the fluid mechanics community: 
'Although breaking and peaking, as well as criteria for 
the occurrence of each, are without doubt contained in 
the equation of the exact potential theory, it is intrigu- 
ing to know what kind of simpler mathematical equation 
could include all these phenomena' [l9(. The short wave 
limit of the CH equation is the Hunter-Saxton (HS) equa- 
tion 



uu x 



0. 



(2) 



obtained from {T]) by taking It describes 

the propagation of waves in a massive director field 



trav- 



flEHEi for u = 0: the peakon u(x,t) = ce~\ x - ct \ 
elling at finite speed c ^ is smooth except at its crest, 
where it is continuous, but has a jump discontinuity in its 
first derivative. The peakons replicate a characteristic of 
the travelling waves of greatest height - exact travelling 
solutions of the governing equations for water waves with 
a peak at their crest [l6|, 17, 18| whose capture by simpler 
approximate shallow water models has eluded researchers 
until recently [191 ]. A further remarkable property of the 
CH equation is the presence of breaking waves (i.e. the 



of a nematic liquid crystal 22|, with the orientation 
of the molecules described by the field of unit vectors 
(cosu(a;,t), sinu(x,t)), where x is the space variable in a 
reference frame moving with the linearized wave velocity, 
and t is a 'slow time variable'. 

The equations ((T|), @ admit many integrable mul- 
ticomponent generalizations [23j, the most popular of 
which is 



wit + 2u x m + um x + 



OPPx = 0, 
Pt + {up) x = 0, 



(3) 
(4) 



where m = a\u — u XXl a = ±1 and o\ = 1 or, in the 
'short wave' limit, cr\ = 0. (The CH equation can be 
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obtained via the obvious reduction p = 0.) This system 
appears originally in 24 1 and its mathematical properties 



have been studied further in many works, e.g. [23l . l25l . 
26 . 27, [2^|. The system is integrable - it can be written 
as a compatibility condition of two linear systems (Lax 
pair) with a spectral parameter (: 

*«»= (-<rCV + Cm+^)*, 
It is bi-Hamiltonian, the first Poisson bracket 



{^1,^2} = 



- — (mo + omj- h - — po— h ——op- — 

dm dm dm dp dp dm 



dx. 



corresponding to the Hamiltonian 
2 



H 



(um + ap 2 )dx 



and the second Poisson bracket 



—(9 -d' i )— + ^d— 
dm dm dp dp 



dx 



corresponding to the Hamiltonian 

1 

2 



Ho 



(oup 2 + u 3 + uu 2 )dx. 



There are two Casimirs: J pdx and J mdx. 

In what follows, we are going to demonstrate how the 
system ©, Q arises in shallow water th eory . We start 
from the Green-Naghdi (GN) equations [H, [2!| , which are 
derived from the Euler's equations under certain assump- 
tions, as follows. Consider the motion of shallow water 
over a flat surface, which is located at z — with respect 
to the usual Cartesian reference frame. We assume that 
the motion is in the x-direction and the physical vari- 
ables do not depend on y. Let h be the mean level of 
water, a - the typical amplitude of the wave and A - the 
typical wavelength of the wave. Let us now introduce the 
dimensionless parameters e = a/h and d = h/X, which 
are supposed to be small in the shallow water regime. 
The variable u(x, t) describes the horizontal velocity of 
the fluid, j](x, t) describes the horizontal deviation of the 
surface from equilibrium, all measured in dimensionless 
units. The GN equations 



U t + £UU X 

TH + [Ml 



<5 2 /3 
l+erj 
er,)]x = 0. 



[(1 + er]) 3 (u xt + euu x 



are obtained under the assumption that at leading order 
u is not a function of z. The leading order expansion 
with respect to the parameters e and d 2 gives the system 



euu x + rj x = 0, 
rjt + [(<it(l + eri)] x = 0. 



(5) 
(G) 



One can demonstrate that the system ([5]), ([6]) can be re- 
lated to the system ([3]), (HJ in the first order with respect 
to e and d 2 . Indeed, let us define 

P = l + \ei 1 -^e 2 (u 2 +1?). 

The expansion of p 2 in the same order of e is 

p 2 = 1 + en - -e 2 u 2 . 

4 

With this definition it is straightforward to write ([5|) in 
the form 

( 52 \ 3 1/ 2x 

u - —u xx + -euu x + -(p )x = u, 
\ 6 / 1 2 e 

or, introducing the variable m = u— ^5 2 u xx , at the same 
order (i.e. neglecting terms of order eS 2 ) 



m t + emu x + \eum x + -(p 2 ) x = 0. (7) 
2 £ 

Next, using the fact that w t « — r\ Xl r\ t ~ —u x , from the 
definition of p we get pt = \^"Qt + \^ 2 ij]u) x . With this 
expression for p t and with p w 1 + eu, equation © can 
be written as 



Pt + ^(pu) x = 0. 



(8) 

The rescaling u -> fu, a; -> ^x, f -> in Q, © 
gives ©, ([!]) with a = <j\ = 1. The case cr = — 1, 
which is often considered, corresponds to the situation 
in which the gravity acceleration points upwards. We 
mention also that the Kaup - Boussinesq system [3(| is 
another integrable system matching the GN equation to 
the same order of the parameters e, 5 0, [3l[ . Notice that 
in the hydrodynamical derivation of (3)- (4) we require 
that u(x,t) — » and p(x,t) — > 1 as |x| — > 00, at any 
instant 

We will now show that for the system ([3]), ([4]) in the 
hydrodynamically relevant case a = a% = 1 wave break- 
ing is the only way that singularities arise in smooth so- 
lutions. The system admits breaking wave solutions as 
well as solutions defined for all times. In particular, we 
will analyze the traveling wave solutions 

The well-posedness (existence, uniqueness, and con- 
tinuous dependence on data) follows by Kato's semi- 
group theory [32| for initial data uq = it(-,0) g H 3 
and po = p(-,0) such that (p — 1) G H 2 If 
T = T(uq,pq) > is the maximal existence time, then 
the integral of motion 



[u 2 + u\ + (p- l) 2 ]dx 



(9) 



ensures that u(-,t) is uniformly bounded (i.e. for all val- 
ues of x £ R and all < t < T) in view of the inequality 

sup \u(t,x)\ 2 < - f(u 2 + u 2 x )dx. (10) 
iei 2 J 
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Considerations analogous to those made in [2jj for a sim- 
ilar system show that the solution blows up in finite time 
(i.e. T < oo) if and only if 



liminf {u x {t,x)} = — oo, 



(11) 



which, in light of the uniform boundedness of u, is inter- 
preted as wave breaking. 

To show that wave breaking occurs, we introduce the 
family {(f(-, i)}te[o,T) °f diffeomorphisms tp(-,t) : R — > K 
defined by 

d t (p(x,t) — u(tp(x,t),t), ip(x,0)=x, (12) 
and we denote 

M(x, t) = u x (cp(x, t),t), y(x, t) = p(ip(x, t), t). 
Consider now initial data satisfying po(0) — and 

^(0)<-2(|| Uo ||2 + ||p-l||2) 1/2 . (13) 

Noticing that (1 — d 2 )^ 1 f = p * f (convolution) with 
p(x) = — e - ' 31 ', and applying the operator (1 — d 2 )^ 1 to 
(3), we get 

u t + uu x + p * (u 2 + i u\ + - p 2 ) = 0. 

Applying now d x and using the identity d 2 p*f = p*f—f, 
we obtain 

u tx +uu xx +^ul = ^(u 2 +p 2 ) -p* (u 2 +^u 2 x +^p 2 ). (14) 
This equation in combination with (|12[) yields 

8 t M(t, x)+^M 2 (t, x) < X - (u 2 (<p(x, t), t)+ 7 2 (a:, *)) • ( 15 ) 

On the other hand, from (|12p and (4), we obtain 

d a = - 7 M. (16) 

Since 7(0,0) = we infer that 7(0, t) = for < t < 
T. The relation $T5fy together with ©, flT0J) ensure that 
Au 2 {ip{x,t),t) < M 2 (0,0). But then fHfy yields 

d t M(0, t)<- ^M 2 (t 7 0) for < t < T. 
As M(0,0) = u' (Q) < 0, this implies 

K(o) 



M(0,t) < 



i + u' o {0)t 



-oo 



in finite time. 

However, not all solutions develop singularities in fi- 
nite time. For example, if the initial data is suffi- 
ciently small, then the solution evolving from it is de- 
fined for all times. More precisely, let a G (0, 1) and 



assume that |1 — po(x)\ < a for all x G R, while 
1 1 1 1 2 + \\po — 1 1 1 o < a - Then the corresponding solu- 
tion to ((3|), J4} is global in time. Indeed, if the maximal 
existence time were T < oo, then for some Xq G R we 
would have liminf t |T M(t, xq) — — oo. We now show 
that this is impossible. First, notice that 



u 2 (x,t) < a 2 , (x,t) G 



[0,T). 



(17) 



by (JH) , (fl7J| . The integral of motion ([9]) also ensures that 
||p(-,*)-l||o < ol on [Q,T). Using this and Q, ®, (flTf. 
we get 



< p * (u 2 + - [u 2 



P 2 ]) 



< \p * " 2 + \ (" 2 + u l + (P- 1 f ) + P * + \p * 1 



; ci 1 



a - 



1 



on Rx [0, T). We now infer from flT2J| , |HJ) that at x = x , 

(18) 



M t = -\M 2 + \ 1 2 -f{t), fe[0,T), 



with the continuous function / : [0, T) — > [0, oo) 
bounded, i.e. there is a constant fco > such that 
< 7(t) < /c on [°, T )- The solution ( 7 (t), M(t)) of 
the nonlinear system (fl7q|) . (fT5|) with initial data 7(0) = 
po(^o) > and M(0) = Uo(xo) supposedly blows up in 
finite time as liminf t ]T M(t) = —00. However, notice 
that (HH) and 7(0) > ensure y(t) > for all t G [0, T). 
Thus we may consider the positive function 

w (t)= 1 (0Mt) + ^[l + M 2 (t)} 

for t G [0,T). Using JTBJ), (JUJ) we get 



w'(t) 



< [k + l}A {1 + M 2 )<[k + l}w(t) 

for all t G (0,T). Thus 

w(i) < io(0) exp([fc + 1]<) 

on [0, T) and this prevents blowup. The obtained con- 
tradiction shows that the solution is defined globally in 
time. 

Global existence is however not confined to small am- 
plitude waves, as we shall see now by establishing the 
existence of traveling waves of large amplitude: solutions 
u(x,t) = ip(x — ct), p(x,t) — £(x — ct) traveling with 
constant wave speed c > 0. To find whether such solu- 
tions exist, notice that with the previous Ansatz equation 
(4) becomes £'(c — ip) = and the asymptotic limits 
i/j(x) — > and £(x) — > 1 as \x\ — > 00 yield 
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Thus (3) becomes a differential equation solely for the 
unknown if). Integrating this equation on (—00, x] and 
taking into account the asymptotic behaviour of tp, we 
get the equation 



cip + ctfj" + -ip 2 



1 r 2 I 

vv 2 W ' 2(c-ip) 2 2 



Multiplication by ip' and another integration on (— 00, x] 
leads to 

(V) 2 -^ 2 ) (c-^) + ^=^ + c, 
recalling the decay of ip far out. Thus 



th 2 



(c-ipy 



(c-V-l)(c-V + l)- (19) 



The asymptotic behaviour ^(2;) — > as — > 00 yields 
now the necessary condition c > 1 for the existence of 
traveling waves. For c > 1, a qualitative analysis of (fT5|) 
shows that < ?/> < c — 1. Thus nontrivial traveling 
waves exist only for c > 1, in which case both ip and £ 
are smooth waves of elevation with a single crest profile 
of maximal amplitude c — 1, respectively c. It is possible 
to find explicit formulas for the traveling waves in terms 
of elliptic functions cf. 34]. Due to the integrability of 



the system, we expect the solitary waves to interact like 
solitons. 

Notice the absence of peakons among the solitary wave 
solutions. However, there are peakon solutions of the 
'short wave limit' equation a\ = 0. Although this limit is 
not covered by the presented hydrodynamical derivation, 
we will describe briefly the construction of the peakon 
solutions, since these are interesting by themselves. The 
limit CTi = is a two component analog of the Hunter- 
Saxton equation. Such system is a particular case of the 
Gurevich-Zybin system [35(, which describes the dynam- 
ics in a model of nondissipative dark matter [3^ |. 

The peakon solutions have the form 

N 

m(x,t) = ^m k (t)5(x - x k (t)), (20) 
fe=i 

1 N 

u(x,t) = --^2m k (t)\x-x k (t)\, (21) 



fc=i 



N 



P(x,t) = ^2p k {t)9(x - x k (t)), 



(22) 



fc=i 



where 8 is the Heaviside unit step function. The asymp- 
totic behaviour p(x, t) — > for x — > 00 and the condition 
/ m dx = (recall that m = —u xx ) lead to 



N 



^2 mi = ^2 pi = 0, 



N 

E 

1=1 



pi = 



in terms of the new complex variable p k = m k J rip k . The 
substitution of the Ansatz (J20J) - (J22J) into ©, (@}, under 
the assumption that x\(t) < 2:2 (i) < . . . < xnU) for all t, 
(a condition holding for the peakons of ([2]) cf. [37]]) gives 
the following dynamical system for the time-dependent 
variables: 



dx k 
dt 

dpk 
dt 



Pk 



1 N 

-Y 

2 ^ 

i=\ 

N 



mi \xi - x k \, 



(23) 
(24) 



with the convention sgn(0) = 0. The integrals for this 
system can be obtained from the integrals of (J5J), (JU) 
(available in [23j) by substituting the expressions (|20| . 
(|21[). ([22l . It is convenient to write the system in terms 
of the new independent variables 



with k = 1,2, 
system is 



x k +i -x k , M k = pi + . . . + p k , 

N — 1. The Hamiltonian of the new 



JV-i 



the equations 
dA fc 



-Re(M fc )A fe , 



dM k 



dt v dt 

are Hamiltonian with respect to the bracket 



{A k ,M l } = 



Hk: 



in which the bar stands for complex conjugation. These 
equations integrate immediately: 



M k (t) 
A fc (t) 



t/2 + c k ' 



A fe (0) 



(t/2 + c M ) 2 +c 2 ka 



- k ,i 



where c k = c ki \ + ic k .2 = ~Mj7 1 (0) is a complex constant 
with real and imaginary parts Cfc,i and c k ^. respectively. 
Notice that the large time asymptotics 



M k 



r 



t l 



1=1 



1=1 



are the same as those for the peakons of the Hunter- 
Saxton equation (2) when pk = (see (37j). 
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